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FOREWORD 


This  final  report  (Unclassified  Title)  "Space 
Encounter  Modeling  Study"  describes  research  conducted 
by  the  IIT  Research  Institute  under  Contract  AFQ8(635)- 
4203,  IITRI  Program  K6064.  The  work  period  was 
28  Feb  1964  to  28  Feb  1965.  The  program  was  administered 
by  the  Air  Force  Armament  Laboratory  (ATWR),  Eglin  Air 
Force  Base,  Florida.  Mr.  R.  L.  Hill,  Jr.  was  the  AFATL 
Project  Officer. 

This  report  comprises  four  volumes: 

I  Model  Formulation  (Unclassified) 

II  Mathematical  Models  (Unclassified) 

III  Sample  Problems  (Secret) 

IV  Computer  Program  (Unclassified) 

The  IITRI  Proiect  team  included  L,  A.  C.  Barbarek 
(Project  Engineer),  M.  Liveright(Statistician  and 
Progranmer) ,  E.  P.  Bergmann  (Warhead  Evaluation)  and 
T.  Martin  (Fire  Control). 

Information  in  this  report  is  embargoed  under  the 
Deparment  of  State  International  Traffic  In  Arms  Regulations, 
This  report  may  be  released  to  foreign  governments  by 
departments  or  agencies  of  the  U.  S,  Government  subject  to 
approval  of  the  Air  Force  Armament  Laboratory  (ATWR) ,  Eglin 
Air  Force  Base,  Florida,  or  higher  authority  within  the 
Department  of  the  Air  Force.  Private  individuals  or  firms 
require  a  Department  of  State  export  license. 


This  report  has  been  reviewed  and  is  approved. 

DAVID  K.  DEAN,  Colonel,  USAF 
Chief,  Weapons  Division 


ABSTRACT 


Studies  conducted  to  provide  the  Air  Force  with 
mathematical  models  of  probable  space  tactical  situa¬ 
tions  are  described  in  this  four-volume  report.  The 
need  for  these  models  was  established  in  order  that  the 
Air  Force  have  means  of  guiding  the  logical  and  timely 
development  and  evaluation  of  nonnuclear  space  warhead 
concepts. 

The  model  described  can  determine: 

•  Warhead  fuzing  and/or  firing  time  (projected 
warheads)  for  a  constant  velocity  intercept 
with  a  given  miss-distance  or  miss-distance 
distribution. 

•  Warhead-target  intercept  error  dispersion  as 
a  function  of  warhead  fuzing  and/or  firing 
time  for  any  combination  of  intercept  param¬ 
eter  errors. 

•  Probability  of  hit  of  a  given  size  vulnerable 
volume  by  a  uniform  random  pattern-type 
warhead. 

•  Cver-all  probability  of  hit  that  is  optimum 
for  a  given  miss-distance  distribution. 

•  Overfall  probability  of  hit  using  a  fuzing 
logic  based  upon  a  given  value  of  some  inter¬ 
cept  parameters  (i.e. ,  line-of-sight  range, 
line-of-sight  bearing  rate) . 
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INTERCEPT  GEOMETRY  MATHEMATICS 

The  first  problem  in  evaluating  a  deli very -warhead 
system  is  to  compute  the  position,  velocities,  and 
angles  at  the  ideal  intercept  point  from  the  initial 
launch  parameters.  The  FLYCUM  subroutine  performs  this 
computation. 

Figure  1  shows  the  angles  and  distances  that  must 
be  considered  in  this  problem.  The  initial  parameters 
are  denoted  by  a  1  subscript,  whereas  the  intercept 
parameters  are  denoted  by  a  2  subscript. 

The  given  parameters  are : 

R^  initial  range  of  the  target 
target  bearings 
E^  target  elevation 
V. ___  target  speed 

R?  y  line-of-sight  to  velocity  angle 
B  horizontal  to  velocity  angle  in  the  (i,p) 
or  reduction  plane 

V  the  speed  of  the  delivery  vehicle,  or 

par  particle  cloud. 

The  output  parameters  that  must  be  computed  are: 

R^  initial  range  rate 

initial  angle  (sweep)  rate 
x lit  flight  to  the  intercept  point 

R2  target's  intercept  range 
B  bearing 

Ej  elevation 

R2a  line-of-sight  to  velocity  angle 
B^  horizontal  to  velocity  angle  in  the 
reduction  plane 
R2  intercept  range  rate 

A2  intercept  angle  rate 

angle  between  the  two  lines  of  sight. 

Using  the  cosine  law,  the  intercept  geometry  can  be 
expressed  by  the  following  equation  which  may  be 
solved  for  the  time-of- f I ight. 
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In  the  situations  where  this  equation  has  two  solutions, 
the  minimum  nonnegative  solution  is  chosen.  If  no  non¬ 
negative  solutions  occur,  then  there  can  be  no  inter¬ 
cept  and  the  problem  is  terminated. 

After  the  flight  time  has  been  established,  it  is 
necessary  to  determined  the  position  of  the  target  at 
intercept.  It  is  useful  to  consider  a  coordinate  sys¬ 
tem  (X,Y,Z)  such  that  the  X  and  Y  axes  lie  in  the 
horizontal  plane  and  the  X  axis  is  along  the  motion 
of  the  initial  line-of-sight  in  the  horizontal  plane. 

In  this  system,  the  target  coordinates  are: 

X  =  R,  *  cos  (E, )  , 

Y  =  0|  and  (2) 

Z  =  Rx  *  sin  (E.^  . 

Relative  to  this  coordinate  system  the  velocity  compo¬ 
nents  of  the  target  are: 


cos  (R^a)  *  cos  (E^) 


targ  *  ^-sin  (Rla)  *  sin  (E^  *  sin  (Bla)j  , 

vtarg  *  sln  (Rla>  *  cos  (Ela*  >  and 

*  I  ~OB  (Rla>  *  sin  (E1>  \ 


targ  l+sin  (Bla)  *  sin  (Rla)  *  cos  (E1) 


It  is  easy  enough  to  determine  the  target  position  after 
a  time  of  Tfiit  and,  thus,  to  determine  the  intercept 
position  in  the  R,  B,  E  system  as  show  below. 


* 

This  symbol  represents  the  operation  of  multiplication. 
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= 


'{X1 +  *1  *  'W  +  (Y1 +  Y1  *  w 

2 


+  (Z1  +  Z1  *  Tflit> 


B„  = 


+  tan 


/  Y  +  Y  *  T 

-1  p  i  Tf  lit' 


\xi  +  *i  *  Tint I 


E„  =  tan 


-1 


Z  +  Z  *  T 
*1  1  flit 


M  +  *i  *  W2  +(Yi +  h  *  Tfu 


Once  the  position  of  the  target  at  intercept  has 
been  determi  ed  it  is  possible  to  compute  the  velocity 
angles  relative  to  this  system.  Since  the  velocity  com 
ponents  are  already  expressed  in  the  (X,  Y,  Z)  systems 
it  is  just  necessary  to  perform  a  transformation  into 
the  (r2,  i2j  Poi  system.  As  Figure  1  shows,  this  is 
done  by  a  rotation  around  the  Z-axis  of  dB  =  B»  -  B, 
followed  by  a  rotation  around  the  new  Y-axis  by^E_. 

This  transformation  is  expressed  in  matrix  notation  by 
the  equations: 


/rA  /  cosfE^  jOjSintEjA 

*2  f  0  >1>  °  * 

\p2]  \-sin(E1)  J0,cos(E1)/ 


-sin(DB) ,cos(DB) ,0 
0  ,  0,1 


\  1 

/  a\ 

x 

★  I 

• 

Y 

1 

The  intercept  velocity  angles  may  then  be  deter¬ 
mined: 


R2a  =  tan 


-1 


(<i2  +  P2)1/2/ft2V  B 


2a  P2/R2 


The  rates  can  easily  be  computed  from  the  data 
already  available. 


R  =  V  *  cos  (R  ) 
n  targ  na 


£  -  -targ  1 

n  R 


sin  (R_„) 


na 


n 


where  n  -  1  or  2. 


(a)  At  initial  target  position 


Figure  1  Angles  and  distances  in  target  flight 
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Figure  1  (Cont.)  Angles  and  distances  in  target  flight 
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Figure  1  (Cont.)  Angles  and  distances  in  target  flight 


Launcher  Location 


Figure  1  (Concluded)  Angles  and  distances  in  target 
flight 


Finally,  the  angle  between  the  two  lines-of-sight 
can  be  domputed  using  the  cosine  formula: 


D 

a 


i  i  *  , _ _  oiui 

tan'1  A  tar;  - 

\  Rl/T£lit  vtarg 


Thus,  all  the  required  intercept  parameters  may  be 
computed. 
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ERROR  DISTRIBUTION  MATHEMATICS 


After  computing  the  intercept  situation,  it  is 
necessary  to  compute  the  effect  of  the  various  errors 
in  launching  the  warhead  on  the  position  of  the  war¬ 
head  centroid  relative  to  the  target  at  the  nominal 
intercept  time.  The  ERRCM  and  INDVAR  subroutines  are 
the  programs  that  generate  the  distribution  of  the 
centroid  position. 

Three  general  classes  of  errors  may  affect  the 
position  of  the  warhead  centroid.  First,  there  may  be 
errors  in  measuring  the  position  or  velocity  of  the 
target  by  the  warhead  launching  system.  These  errors 
cause  the  warhead  to  be  launched  along  a  path  different 
from  the  ideal  path,  and  thus  contribute  to  the  possi¬ 
bility  of  missing  the  target.  Second,  there  may  be 
errors  in  launching  the  warhead  along  its  flight  path. 
These  also  contribute  to  missing  the  target.  Third, 
there  may  be  errors  in  the  position  of  a  specific  par¬ 
ticle  relative  to  the  flight  path  of  the  delivery 
vehicle. 

Any  single  error  (if  it  is  reasonably  small  com¬ 
pared  with  the  intercept  parameters)  will  cause  an 
error  in  the  position  of  the  centroid  of  the  particles 
of  the  warhead.  This  error  will  lie  along  a  single 
line.  For  example,  an  error  in  the  launch  bearing  will 
cause  an  error  in  the  centroid  that  is  horizontal  and 
is  perpendicular  to  the  line-of -sight  Lo  the  target  at 
nominal  intercept  time.  If  the  error  is  normally  dis¬ 
tributed,  then  its  effect  will  be  a  normal  distribution 
of  centroid  positions  lying  along  this  same  line.  The 
other  errors  lead  to  the  same  type  of  linear  normal 
distributions.  If  the  various  errors  are  all  independent, 
then  they  may  be  combined  by  expressing  them  all  in  a 
single  coordinate  system  and  adding  their  respective 
biases,  variances  and  covariances.  The  results  of  this 
combination  is  a  three-dimensional  normal  distribution. 

The  errors  in  measuring  the  target  posit > an  may 
arise  from  a  number  of  sources.  These  errors  lie  along 
the  range,  bearing  elevation  (R,  B,  E)  system  or  the 
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line-of-sight  (r,  i,  p)  system  at  the  time  when  the 
warhead  is  launched  (Figure  2) .  In  the  case  of  these 
errors  in  measuring  the  target  position  the  line-of- 
sight  system  is  the  reduction  plane  system. 

A  target  error  may  have  one  of  three  different 
effects,  and  errors  must  be  provided  for  each  of  these 
effects.  First,  it  may  affect  the  measurement  of  range, 
bearing,  or  elevation,  or  the  measurement  of  the  rate 
of  any  of  these  variables.  These  errors  are  denoted  by 
the  following  symbols: 

Measured  range  error  ^  he 

Measured  range  rate  error  = 

Measured  bearing  error  =  ®le 

Measured  bearing  rate  error  ~  B-^e 

Measured  elevation  error  =  E^e 

Measured  elevation  rate  error  =  E. 

ie 

Second,  the  error  may  affect  the  position  measurement 
in  such  a  way  that  an  error  occurs  in  a  random  direction 
of  the  (i,p)  plane.  This  will  cause  an  angular  error 
in  the  (i,p)  plane  denoted  by: 

Measured  angular  error  =  Xle* 

Third,  the  error  may  cause  an  error  in  the  i  or  p 
direction  that  is  measured  in  linear  units,  e.g. ,  feet. 
These  types  of  errors  are  denoted  by  the  symbols: 

Measured  i  linear  error 
Measured  i  linear  rate  error 
Measured  p  linear  error 
Measured  p  linear  rate  error 

For  this  analysis,  we  assume  that  these  errors  are 
all  normally  distributed,  and  independent.  The  bias  of 
any  error  A  is  represented  by  B(A)  and  its  variance 
by  V(  A) . 

The  following  equations  show  the  combined  effect 
of  all  of  the  target  errors  mentioned  previously. 
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Launcher 

Location 


Elevation 


Y 


Figure  2  Relative  position  of  target 


B(ple)  +  R  *  B(Ele)  +  Tinfc  *  (B(ple) 
+  R  *  B(Ele)) 


v(rl>  = 


V(Rle} 


V(Rle> 


B(i1)  =  V{ile)  =  (R  *  cos (E) ) 2  *  V(Ble) 

+  Tint2  *  (V(£lei  +  (R  *  COS(E))2  *  v(Bie)) 
+  R2  *  V(\le) 


v(pi) 


V(ple)  *  R2  *  V(Ele)  +  Tlnt2  *  (V(ple) 
+  R2  *  V(Ele))  +  R2  *  V(\le). 


The  range  R  and  elevation  E  in  these  equations 
are  evaluated  at  the  aiming  time  or  time  of  fire,  and 
Tint  time  between  firing  the  warhead  and  its 

intercept  with  the  target. 

The  errors  in  launching  the  warhead  are  similar  to 
those  in  measuring  the  target  position  except  that  they 
will  arise  from  different  causes,  and  they  will  occur 
when  the  warhead  ideally  hits  the  target.  These  errors 
are  denoted  by  the  symbols: 

Launch  velocity  error  =  Vce 

Launch  bearing  error  =  B^g 

Launch  bearing  rate  error  = 

Launch  elevation  error  =  E2e 

Launch  elevation  rate  error  =  £2e* 

As  with  the  previous  error  type,  an  error  may  affect  the 
launch  in  a  random  direction  leading  to  an  angular  error 
in  the  (i,p)  plane  denoted  by: 

Launch  angular  error  =  \  • 
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Finally,  we  may  be  interested  in  a  particle  of  the 
warhead  that  differs  from  the  centroid  of  the  warhead 
in  that  it  has  additional  errors  denoted  as 

Particle  r  velocity  error  =  r 

Particle  i  velocity  error  = 

Particle  p  Velocity  error  =  jr  ,  and 

Particle  angular  error  = 

It  is  possible  to  combine  all  the  errors  that  are 
related  to  the  intercept  coordinate  system  as  those 
discussed  previously  were  combined.  The  following 
equations  show  the  combined  effect  of  all  the  errors  in 
luanching  the  warhead,  and  in  a  particle  position  rela¬ 
tive  to  the  warhead. 


B(r2)  = 


(B(Vce) 


B(R  )  ) 
pe 


B(i2)  =  R  +  cos(E)  *  B (B2e}  +  Tint  *  (B(i  ) 


+  R  *  cos (E)  *  B(B,  ) ) 

4  w 


B(p2)  -  R  +  B(E2e)  +  Tint  *  (B(p  >  +  R  *  *ri2e>) 


V(r2) 


v(i2) 


Tint  *  <V<Vce>  +  V(V)> 

R2  *  <v<\,e>  +  v (Xpg)  +  cos (E) 2  *  V(B2e)) 
+  Tint2  *  <v(ipe>  +  r2  cos<e>2  *  V(B2e)) 


V(p2)  =  R  *  <VUce)  +  v(Xpe)  +  V(E2e)) 
+  Tint2  *  V(Ppe>  +  r2  *  v®2e»> 


Here,  R  and  E  are  evaluated  at  the  intercept  point. 

At  this  point  there  are  two  sets  of  errors,  one 
expressed  in  terns  of  the  launch  system  and  the  other  in 
terms  of  the  intercept  system.  The  next  problem  is  to 
express  one  (the  intercept  system  errors)  in  terms  of 
the  other  (the  launch  system)  so  that  the  two  sets  of 
errors  may  be  combined  into  one  distribution. 
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The  transformation  between  the  intercept  system 
and  the  launch  system  is  the  product  of  three  rotations. 

A  =  First,  a  rotation  around  the  ^-axis  by  an 
angle  of  -E2, 

B  =  Second,  a  rotation  around  the  new  p2*-axis  by 
an  angle  of  -dB  =  - (B2-B2) , 

C  =  Third,  a  rotation  around  the  new  i2"  ="  il  axis 
by  an  angle  of  E^. 


Rotation  C  Rotation  B 


rl 

costEj) ,0,sin(E^) 

cos (dB) , 

-sin(dB) ,0 

il 

= 

0  ,1,  0 

x 

sin(dB) , 

cos (dB) ,0 

Po 

-sintE-^)  jO^osCE^^) 

0  , 

0  ,1 

^  _ 

_  — 

-I 

Rotation  A 


cos(E2) ,0, 

-sin(E2) 

r2 

X 

0  ,1, 

0 

X 

i2 

sin(E2)  ,0, 

cos(E2) 

P2 

_  _ 

It  is  possible  to  combine  the  three  rotations  into  one 
transformation  F,  by  performing  the  multiplication  of 
the  three  matrices.  This  leads  to  the  following  trans¬ 
formation  between  (r2,  i2>  P2^  and  (*1,  i 2i  P2^ • 


Transformation  F 


rl 

= 

frr’fri’frp 

fir’fii’fip 

X 

“  — 

r2 

A2 

h 

fPr ,fpi ,fpp 

P2 

— .  — 

These  coefficients  permit  the  total  bias,  TB (a^) ,  total 
variance,  Tv (a.) ,  and  total  covariance,  Tcv(a^,bi)  to  be 
computed  for  a  and  b  equal  to  r,  i,  or  p.  Thus,  the 
distribution  coefficients  along  each  of  the  launch  axes 
may  be  computed  from  the  equation: 


TB(ai) 


B(al>  +  £ar  *  B(r2>  +  £ai  *  B(12> 
+  £ap  *  B(p2> 
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TV(a1)  =  V(a;L)  4  far2  *  V(r2)  +  fai2  *  V(i2) 

+  £ap2  *  V(p2) 

CV  U1>bl)  =  far  *  £br  *  V(r2)  +  fal  *  £bl  *  V(i2) 

+  £ap  *  £bp  *  V(f2> • 

The  errors  in  the  redirection  will  not  be  consid¬ 
ered  in  the  cloud  models  used.  These  errors  only  affect 
the  time  at  which  a  particle  hits  the  target,  and  only 
indirectly  affect  the  hit  probability.  This  indirect 
effect  is  due  to  the  expansion  of  the  warhead  and,  since 
it  is  assumed  that  the  expansion  speed  is  small  relative 
to  the  relative  speed  of  the  warhead,  the  errors  in  the 
r^-direction  will  be  neglected. 

Thus,  the  distribution  computed  is  projected  on  the 
reduction  (ii,Pi)  plane  by  neglecting  any  of  the  terms 
with  r^  in  them.  Once  the  distribution  is  reduced  to  a 
two-dimensional  normal,  it  is  necessary  to  find  its 
major  and  minor  axes  so  that  the  distribution  may  be 
expressed  as  the  product  of  two  independent  normal  dis¬ 
tributions.  Before  being  rotated,  the  distribution  will 
have  the  form: 

prob  (i,p)  =  K  exp  -  (i2/V,  +  2  1  p  +  p2/V  ) 

1  1  ip  p 

where 

-  Tv(i1) ;  Vp  =  Tv (p^) ;  and  CVip  =  Tcv(i1,p1). 

One  of  the  axes  of  the  independent  distribution  is  along 
the  angle  0  where 

6  =  i  tan"1  (2*CV  /(V  -V. )) 

2  ip  p  i 

The  variance  along  this  axis  is: 

V  .  =  (v  *  cos (0))  4  (v.  *  sin(0)) 
maj  p  i 

+  (2*C . ip*cos (0)  *  sin(0)) 
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and  along  the  other  axis  is: 


Vmin  =  (Vp  *  +  (vi  *  cos(0) 

-  (2  *  CVip  *  cos(e)  *  sin (0 ) ) . 

If  Vma-j  is  greater  than  Vm^n,  then  0  is  left  alone. 
Otherwise ,  90  deg  is  subtracted  to  generate  a  new  0 ,  and 
Va  and  are  interchanged.  This  yields  the  biases 

along  the  major  and  minor  axes: 


Bmaj  “  ra^i1)  *  cos(0)  +  TB(P;l)  *  sin(0) 


Bmin  =  *  sin(0)  +  TB(P;l)  *  cos(0) 


The  five  parameters  0,  Vmaj  ,  Bmaj  ,  Bmin 

fully  describe  the  distribution  of  the  intercept  between 
the  warhead  centroid  and  the  reduction  plane.  Any  fur¬ 
ther  analysis  of  the  hit  probabilities  depends  on  the 
warhead  and  target  characteristics. 
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THREE 


PROBABILITY, OF  HIT  MATHEMATICS 


1 .  Uniform  Random  Pattern 


After  the  flight  and  distribution  parameters  of  the 
warhead  centroid  have  been  computed,  it  is  possible  to 
compute  the  probability  of  hitting  a  relatively  small 
target  with  an  expanding  disk  of  particles.  The  PROCM 
routine  generates  the  probability  of  hit  for  the  U.R. 
warhead. 

The  first  problem  in  this  analysis  is  to  determine 
the  radius  of  the  disk  of  particle  as  it  passes  through 
the  reduction  plane.  The  radius  of  the  particle  disk 
will  be  determined  by  its  expansion  speed  V  and 
its  time  of  flight  T..  :  exp 

R  =  T.  .  *  V 
ave  int  exp 

Since  the  disk  will  usually  not  pass  through  the 
reduction  plane  head  on,  it  will  leave  an  elliptical 
trace.  If  the  angle  between  its  velocity  vector  and 
the  perpendicular  to  the  reduction  plane  is  Dtrg  the 
disk  will  leave  an  elliptical  trace  having  a  major 
axis  Rave  and  a  minor  axis  (Rave)  (cos  (Dtrg)). 

The  program  supplied  by  Dahlgren  (Ref.  1)  permits 
the  calculation  of  the  probability  that  the  center  of 
the  vulnerable  area  is  covered  by  the  particle  disk. 

To  use  this  program  a  circular  particle  disk  must  be 
used,  and  the  error  distribution  must  be  expressed  as 
an  independent  bivariate  normal  distribution.  Thi3 
may  be  achieved  by  transforming  the  coordinater  system 
of  the  reduction  plane.  Figure  3  shows  the  relevant 
angles  and  parameters  in  the  reduction  plane. 

The  first  step  in  this  transformation  is  to  express 
the  normal  distribution  in  terms  of  the  disk  coordinates. 
This  requires  a  rotation  of  the  coordinate  system  of  the 
aiming  dj  stri'bution  by  an  angle  of  rot  where  rot  -  (Hor 
to  aiming  distribution  major  axis  -  Hor  to  distribution 
minor  axis).  As  discussed  in  Cramer  (Ref.  2)  any  linear 
transformation  between  X,  Y  coordinates  and  X? ,  Y1 
coordinates  may  be  expressed  by: 
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i.e. ,  X' 

Y' 


ax  +  by 
cx  +  dy 


For  such  a  transformation, 


bias (X1 ) 
Ibias  (Y1 ) 


M  * 


bias (X) 
bias (Y) 


Particle 

Disk 


Figure  3  Reduction  Plane 


r 


and 

var(X')  ,  covar {X'Y'J 
covar (X'Y'),  var(Y’) 


M* 


var(X)  , 
covar (XY) , 


covar (XY) 
var (Y) 


T 


L8 


but,  for 


a  ,  b 

M  =  . 

c  ,  d 


var(X)  covar(XY)  ^ 

M*  covar(XY) ,  var(Y)  *M' 

[a2  *var (X)  +  b2  *  var(Y)  ac  *  var(X)+  bd  *  var(Y)l 


+  2ab  *  covar(XY) 


+  (ab  +  be)  *  covar(XY) 


1C)  +  bd  *  var(Y)  c 
bd)  *  covar(XY) , 


2  *  var(X)  +  d2  *  var(Y) 


*  var(X)  +  bd  *  var(Y)  c  *  var(X)  +  a  * 
+  (ad  +  bd)  *  covar(XY) ,  +  2cd  *  covar(XY) 


for  a  rotation  of  <  rot 


a  =  cos (rot) 
d  =  cos (rot) 


b  =  +sin(rot) 


c  =  -sin(rot) 


Thus,  in  rotating  to  the  disk  system 

Bias (A)  =  B  .  *  cos (rot)  +  B  *  sin (rot) 

Maj  Min 

Bias(B)  =  BM  ,  *  sin  (rot)  +  B...  *  cos  (rot) 

Maj  Mm 

Var(A)  =  V .  *  cos2  (rot)  +  V  *  sin2  (rot) 
Maj  2  Mm  o 

Var(B)  =  V  aj  *  sin*  (rot)  +  V*  *  cos  (rot) 
M  J  Min 


Covar(A,B)  =  cos (rot)  *  sin(rot)  *  (VMJ  -  V  ) 

ru.il  Maj 


After  the  distribution  has  been  expressed  in  the 
disk's  system  it  is  possible  to  "elongate"  ftie  A  coor¬ 
dinate  so  that  the  disk  may  be  considered  to  have  a 
circular  trace  in  the  reduction  plane.  This  is  done  by 
defining  a  new  coordinate  system  A'B'  so  that: 
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A*  =  - rfe: — 7  *  A  and  B'  =  B, 

cos (Dtrg)  ’ 

using  the  previous  notation,  this  transformation  is 
obtained  by  setting 


1 

cos (Dtrg) ’ 


b  =  c  =  0,  and 


d  =  1. 


and  the  distribution  hao  new  parameters  as  follows 

Bias  (A1)  =  - — r  Bias  (A) 

cos (Dtrg) 

Bias (B 1 )  =  Bias (B) 

Var(A')  =  - -  Var(A) 

cos^ (Dtrg) 

Var(B')  =  Var(B) 

Covar(A'B' )  =  — (4'rg7  CovarlA.B) 

Next  it  is  necessary  to  find  the  principal  axes  of 
this,  distribution,  and  express  all  the  parameters  of  the 
distribution  in  terms  of  their  coordinate  system.  Since 
this  is  the  same  transformation  required  in  determining 
the  major  and  minor  axes  of  the  original  distribution, 
no  further  discussion  is  needed  of  this  step. 

At  this  point  the  data  are  in  a  form  such  that  the 
Dahlgren  (Ref.  1)  program  may  be  used  to  compute  the 
probability  that  the  center  of  the  target  is  within  the 
particle  disk. 

This  program  computes  the  probability  that  a  disk 
of  radius  Raye  wij-1  cover  a  point  if  the  center  of 
the  disk  is  distributed  as  a  bivariate  normal.  This  is 
given  by  the  integral; 
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'R 


ave 


P (Cover)  =  j  dx 

i 


^-R 


ave 


where 


exp 


lx-x)2,  (Y-Y)2 
2  2 
.  2ox  2ay2 


dy 


a  and  cr  are  the  standard  deviation  of  the 
^  distribution,  and 

x  and  y  are  its  biases. 


This  computation  results  in  the  probability  that 
the  center  of  the  target  is  covered  by  the  particle 
disk. 


Since  the  vulnerable  region  is  assumed  to  be  rela¬ 
tively  small  compared  to  the  size  of  the  particle  cloud, 
the  probability  that  the  area  is  partially  covered  by 
the  cloud  is  low  relative  to  the  probability  of  being 
completely  covered  or  missed  as  discussed  previously. 

In  this  case,  it  is  possible  to  assume  that  the  area 
is  cohered  when  the  center  is  covered.  When  the  area 
is  covered  by  the  particle  cloud,  it  is  possible  that 
all  the  area  is  not  hit  if  the  particles  of  the  cloud 
miss  the  area.  When  there  are  a  large  number  of  par¬ 
ticles  and  the  area  is  much  smaller  than  the  cloud,  the 
probability  that  K  particles  hit  an  area  within  the 
cloud  may  be  approximated  by  a  Poisson  distribution: 

P (K  hits)  =  ~r  *  e~A- 

In  • 


where 


X  =  q  *  AA/disk  area 
r)  =  number  of  particles  in  disk 
AA  =  area  of  vulnerable  region. 

It  should  be  noted  that  the  area  of  the  disk  is  not  the 
area  of  the  original  circular  disk,  but  of  the  elliptical 
trace  in  the  reduction  plane.  The  probability  of  no 
hits  is  e“^,  therefore,  the  probability  of  one  or  more 
hits  is  1  -  e-^. 


\ 
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It  is  then  possible  to  compute  the  probability 
that  an  area  is  hit  by  multiplying  the  probability 
that  it  is  covered  by  the  cloud  by  the  probability 
that  if  it  is  covered  by  the  cloud  it  will  be  hit. 

P(hit)  =  P (cover)  *  (1  -  e~X)  . 

At  this  point  the  required  value  has  been  computed, 
and  the  computer  program  is  finished  with  the  specific 
intercept  situation,  and  either  considers  another  inter¬ 
cept  or  stops. 

2 .  Continuous  Rod  Mathematics 


To  determine  some  of  the  problems  that  will  be  en¬ 
countered  in  studying  complex  warheads  against  complex 
targets  it  was  decided  to  study  the  problem  of  the  con¬ 
tinuous  rod  warhead  against  a  number  of  vulnerable  re¬ 
gions.  The  warhead  is  an  expanding  solid  ring  that  is 
assumed  to  disable  the  target  if  the  ring  hits  any  of 
the  vulnerable  regions  of  the  target. 

Mathematical  Statement  of  Problem 


In  completely  general  terms,  the  probability  of 
any  event  E  that  depends  on  other  variables  x1?  x2,. 
....,xn  is  defined  by  the  equation: 


Pr  (E) 


xi :x2  ? • • *xn 
such  that  E 

occurs 


p(xl,x2, 


,xx)  dxl?dx2 


.  dx 


n 


An  equivalent  way  of  expressing  this  is  to  use  the  indi¬ 
cator  function 


b  (***)  = 


1  when  ***  is  satisfied 
0  when  ***  is  not  satisfied 


Then,  it  is  possible  to  express  the  probability  of  an 
event  E: 


Pr<E>  =aU  f  '>(E(x1,x2,.,xn)occurs)P<xi>x2’--xn)dxldx2-dx„- 


X1 ’x2 ’ *  *  * ,xn 
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In  evaluating  a  warhead,  the  event  E  to  be  con¬ 
sidered  is  that  of  intercepting  the  target  with  the 
warhead.  The  event  of  not  intercepting,  -E ,  may  also 
be  useful  in  further  work.  The  four  random  variables 
that  affect  the  event  are  i  and  p,  the  position  of 
intercept  of  the  warhead  center  with  the  reduction  plane, 
Ve,  the  velocity  of  the  expansion  of  the  warhead;  and 
T^,  the  time  of  detonation  of  trie  warhead  before  it 
intercepts  the  reduction  plane.  These  last  two  vari¬ 
ables,  Ve  and  T^,  may  be  combined  to  give  the  size 
of  the  ring,  Sr,  at  intercept: 

Sr  -  <Ve>  <  <Ta) 

=  (V  )  *  (Vw/Dd) 

G 

where  Vw  is  the  velocity  of  the  warhead  relative  to  the 
target  and  Dd  is  the  distance  from  warhead  to  reduction 
plane  at  detonation.  Then: 


Pr (hit) 


f  b (hit/i ,p,s) 
i,p,s 


p(i,p,s)  di 


dp  ds 


(1) 


The  p(i,p,s)  term  depends  on  the  distribution  of  errors 
in  aiming  and  detonating  the  warhead.  If  the  distribu¬ 
tion  of  ring  radii  is  independent  of  the  center  location 
and  p(i,p)  is  a  bivariate  normal,  then,  by  proper  def¬ 
inition  of  the  (i,p)  direction,  the  distribution  may  be 
broken  into  three  independent  terms: 

-p(i,p,s)  =  p  (i)  *  p(p)  *  p(s).  (2) 

The  function  6(hit/i,p,s)  is  somewhat  more  difficult 
to  evaluate.  This  term  depends  on  the  shape  of  the  ring 
intersection  with  the  reduction  plane  and  on  the  vulner¬ 
able  regions.  The  actual  shape  of  the  ring  is  an  annulus 
of  relatively  small  width,  but  it  may  not  be  perpendicu¬ 
lar  to  the  relative  flight  path  of  the  warhead,  thus, 
may  intersect  the  reduction  plane  in  an  ellipse. 

Consider  the  situation  shown  in  Figure  4.  If  the 
ring  is  at  an  angle  a  with  the  reduction  plane  and 
the  radius  is  s,  then  the  intersection  of  the  ring 
with  the  (X , Y)  plane  is  an  ellipse  given  by 

x  /cos (a)  +  y  =  s  . 
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Figure  4  Intercept  of  ring  and  reduction  plane 


The  indicator  function  for  a  set  of  vulnerable  areas 
is  the  Boolean  sum  of  indicator  functions  for  each  vul¬ 
nerable  region.  This  function  is  one  if  any  vulnerable 
region  is  intercepted  by  the  ring,  and  zero  otherwise. 

Let  us  consider  a  relatively  simple  case  of  one  vulner¬ 
able  region  where  the  ring  intersection  and  the  vulner¬ 
able  area  are  both  circles.  In  this  case  the  volume 
where  the  interception  occurs  is  a  truncated  hollow  cone. 
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If  the  vulnerable  region  is  centered  at  the  point 
(0,0)  with  radius  R,  interception  occurs  if  i,p,  and 
s  satisfy  both  the  equations  (see  Figure  5) . 

i2  +  p2  ;>  s2  -  R2  and  i2  +  p2  ^  s2  +  R2 


Figure  5  Intercept  of  circular  ring 

and  circular  vulnerable  region 


When  more  than  one  vulnerable  region  is  involved, 
the  intercept  occurs  when  the  i,  p,  and  s  values 
satisfy  one  or  more  sets  of  equations  of  the  type  shown. 

Two  factors  complicate  this  analysis  in  the  real 
problem:  (1)  First,  the  ring  intersection  is  an  el¬ 

lipse  rather  than  a  circle,  and  (2)  more  than  one  vul¬ 
nerable  area  is  involved.  The  first  problem  means  that 
instead  of  being  cones  with  circular  cross  sections, 
the  boundaries  of  the  intersection  volumes  have  approx¬ 
imately  elliptical  cross  sections.  The  second,  and 
really  more  difficult  problem,  is  that  a  number  of 
these  cone-shaped  areas  must  be  considered.  The  volumes 
that  must  be  integrated  over  are  thus  relatively  complex, 
and  would  be  difficult  to  characterize,  even  in  the  sim¬ 
ple  case  where  the  ring  intersection  is  a  circle. 
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A  Numerical  Technique  for  Evaluating 
a  Continuous  Rod 


In  the  case  of  where  the  size  s  is  independently 
distributed  relative  to  the  intersection  point  i,p, 
t'e  probability  of  hit  by  a  warhead  may  be  given,  in 
general,  by  combining  evaluations  (1)  and  (2)  to  obtain 
the  equation: 


Pr (hit)  =  j 

s 


1  p 


&(hit/i,p,s)  p(i,p)  di  dp)j  p  (s)  ds 


To  analyze  this  by  numerical  methods  it  is  neces¬ 
sary  to  approximate  the  integrals  by  sums: 


s^s^k^As 
Pr(hit)  =  ^ 

s=so~k1As 


io+k2Ai  po+k3Ap 

^>~  b  (hit/i  ,p,  s)  p  (i  ,p)  AiAp 
\i“^c-k2Ai  p=po-k3Ap  1 


*  p  (s)  As 

The  indicator  function  &(s,i,p)  specifies  which 
combinations  of  s,i,  and  p  result  in  a  hit.  This 
function  is  independent  of  the  distribution  of  these 
parameters.  It  was  decided  to  evaluate  this  indicator 
function  first  since  it  may  then  be  used  to  evaluate  a 
number  of  different  aim-point,  and  ring-size  distribu¬ 
tions.  The  program  will  determine  the  indicator  func¬ 
tion  relatively  early,  and  later  the  program  will  sum 
this  function  over  various  distributions. 

This  program  may  be  divided  into  major  processing 
stages : 

I.  Computation  of  &(s,i,p)  for  a  set  of 
grid  points  in  the  (i,p)  plane. 

II.  Computation  of  Pr (hit/s)  = 

j-  ^  6(s,i,p)  p(i,p)AiAp 

for  any  specific  distribution 
of  grid  points. 
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6  (s,i,p) 
Computation 


Initial  Shape 
Input 


Transformation 
of  Ellipse 
to  Circle 
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of  Test  Points 
(Grid  Points) 


Computation 
of  Distances 
between  Grid 
Points  and 
Vulnerable  Areas 


Cycle  through 
each  grid  point 
and  each  vul¬ 
nerable  area 


P  (hit/s) 
Computation 


Input  of  Centroid 
Distribution  and 
Computation  of  Grid 
Square  Probability 


III  Pr  (hit 


Computation 


©  1 


Computation  of 
Probability  of  Hit 
for  Each  Ring  Radius 


Cycle 
through 
each  grid 
point 


Input  of  Radius 
Distribution  and  Com¬ 
putation  of  Probability 
of  Hit  for  the 
Specific  Distribution 

Return  to  ®@®  or  terminate 
— »  depending  on  tKe  input 


Figure  6  Flow  chart  of  continuous  rod  warhead  program 
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III.  Computation  of  Pr(hit)  for  any  specific 
distribution  of  ring  rsizes. 

These  three  major  processing  stages  are  then 
divided  into  a  total  of  seven  stages.  A  flowchart  of 
these  stages  is  given  in  Figure  6. 

Stage  A:  Initial  Input  Data  (Shape  Data) .  Some  input 
data  are  required  before  any  processing  can  be  done. 

These  data  are  a  description  of  the  vulr  rable  areas, 
the  shape  of  the  ring  warhead  intersection  with  the 
reduction  plane,  and  estimation  as  to  the  size  of  the 
various  distributions,  and  the  accuracy  desired  in  the 
final  answer. 

The  description  of  the  vulnerable  areas  could  be 
given  in  a  number  of  ways.  To  simplify  some  of  the 
computations,  it  was  decided  to  describe  the  vulnerable 
areas  as  a  limited  number  of  polygons  each  with  a 
limited  number  of  sides.  These  data  are  stored  in  the 
following  way: 

iar  =  number  of  areas 

ico(i)  =  number  of  corners  or  sides  cf  each  area 
cx(i , j ) ,cy (i , j )  =  coordinates  of  the  jfch  corner 

of  the  ith  shape. 

The  shape  of  the  ring  intersection  is  given  by  its 
orientation  and  eccentricity.  These  will  be  given  by 
an  angle  and  a  ratic  of  the  major  to  minor  axis: 

rt  =  angle  of  major  axis  with  X  axis 
re  =  major  to  minor  axis  ratio. 

The  size  of  the  aim  points  distribution  and  the 
rings  size  distribution  must  be  estimated  to  determine 
the  grid  size  that  must  be  used  in  later  stages. 

sc  =  The  average  standard  deviation  of  the 
centroid  distribution 
_  ,  /2  (maximum  standard  deviation  1 
'  minimum  standard  deviation} 
sy  =  coordinates  of  the  centroid  aim  point 
sr  =  the  standard  deviation  of  the  ring  radius 
r  zero  =  average  ring  radius. 
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These  values  are  estimated  so  that  the  proper  grid  size 
is  used  in  each  case.  Various  actual  distributions  will 
be  input  later.  If  these  values  are  properly  estimated, 
then  one  calculation  will  suffice  for  a  number  of  dif¬ 
ferent  distributions.  The  accuracy  of  a  specific  cal¬ 
culation  will  depend  on  the  grid  size.  To  reduce  the 
computation  time,  it  is  desirable  to  be  able  to  control 
this  size.  This  is  done  by  specifying  the  number  of: 
divisions  of  the  grid  in  each  direction  in  the  reduc¬ 
tion  plane: 

ngd  =  number  of  grid  divisions  in  each  direction. 

Stage  B:  Ellipse  to  Circle  Transformation.  Since  it 
is  easier  to  deal  with  a  circular  ring,  ths  first  trans¬ 
lation  is  to  distort  the  X,Y  plane  into  an  (i,p)  plane 
in  which  the  ring  intersection  is  a  circle.  This  is 
represented  by  a  rotation  by  -0,  A,  a  stretch  of  Y 
axis  by  re,  B,  and  a  rotation  by  +0,  C. 


_ 

C 

_  E  n 

| 

A 

r 

r 

i 

cos0 ,  sin0 

1  0 

cos0  -sin© 

X 

P 

-sin6.  cos0 

• 

0  e 

• 

[ 

sin©  cos© 

• 

Y 

0  =  rr  and  e  =  re 


This  transformation  should  be  applied  to  the  co¬ 
ordinates  of  the  corners  as  the  vulnerable  areas  and 
the  bias  of  the  centroid  distribution  to  obtain  the  new 
values  in  the  (i,p)  plane.  To  save  storage  the  trans¬ 
formed  values  will  be  stored  in  the  original  locations, 
and  be  denoted  by  the  same  symbols.  This  transformation 
should  be  saved,  since  other  values  must  be  expressed 
in  the  (i,p)  system. 


Stage  C:  Grid  Generation.  It  is  necessary  to  generate 
a  set  of  grid  points ,  each  to  be  used  as  a  location  for 
the  center  of  the  ring.  The  separation  of  these  points 
will  depend  on  the  number  of  grid  divisions  (ngd)  de¬ 
sired,  and  or.  the  average  standard  deviation  of  the  aim 
point  sc.  One  technique  for  choosing  the  separation 
is  to  let  the  probability  of  a  tail  of  the  normal  be 
equal  to  an  average  probability  in  any  specific  slot. 
The  average  probability  of  being  in  any  slot  is  approx¬ 
imately  1/ngd.  A  rough  estimate  of  the  probability 
of  being  in  a  tail  is 


~=l  j  exp  (-x2/2) 
v  2tt 

where  x  =  cutoff/0.  A  relatively  simple  table  can  be 
stored  for  the  cutoff  point,  cut(ngd),  as  a  function  of 
ngd.  Then  the  separation,  sep,  is  given  by 

sep  =  sc  .  c^t (ngd) /ngd. 

The  grid  points  will  then  be  selected  by: 

x  =  sx  ±  k.  *  sep 
gp  1 

Ygp  =  sx  ±  k2  *  sep 

with: 

0  £  k  £  ngd/2 
0  <  k2  <  ngd/2. 

Distance  Determination.  For  each  grid  point,  a  set  of 
radii  that  the  ring  can  have  and  intercept  a  vulnerable 
area  will  be  computed.  Let  us  consider  one  area  with 
corners  (in  the  (i,p)  plane),  cx . ,  cy . ,  and  a  grid 
point  at  g  =  (gx,  gy) .  The  distance1 to  a  corner  is 
given  by: 

2  2  ^//2 
disi  =  ((cxi-gx)  +  (cy^gy)  ) 

The  maximum  of  these  distances  is  the  maximum  radius 
the  ring  may  have  and  still  intercept  the  area.  This 
will  be  called  rmax.  Similarly,  a  minimum  distance 
corner,  tentatively  called  rmj_n,  exists.  Two  further 
considerations  must  be  considered  in  determining  rm^n. 
First,  it  is  possible  that  the  grid  point  is  inside 
the  vulnerable  area.  In  this  case,  rm^n  =  0.  Second, 
the  ring  may  pass  through  a  side  rather  than  the 
nearest  corner  and,  in  this  case,  rmin  is  less  than  that 
based  on  the  corners. 

To  test  whether  the  point  is  inside,  we  use  the 
Jortan  curve  theorem.  Consider  the  infinite  half-line 
drawn  from  the  grid  point  (gx,  gy)  to  the  "infinite" 
point  (gx,  00).  A  line  between  two  corners  (cxj_,  cy^) 
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and  (cxj+i,  cYi+i)  is  crossed  by  the  infinite  half¬ 
line  if,  and  only  if, 

(cx^-gx)  *  (cxi<rl-gx)  <  0 


and 


ay 


(cyi+1-cYi)  * 


CX2  -  gx 
CX2  -  cx^ 


A  point  is  inside  the  area  if  the  infinite  half-line 
crosses  such  boundaries  an  odd  number  of  times. 


If  a  corner  of  the  area  is  on  the  half-line,  then 
an  intersection  of  the  half -line  with  the  boundary 
occurs  when  the  y  coordinate  of  the  corner  is  greater 
than  gy  and  the  two  adjacent  corners  are  not  on  the 
same  side  of  the  half-lines.  Thus,  it  is  possible  to 
determine  when  a  grid  point  is  inside  an  area. 

The  second  consideration  occurs  when  the  grid 
point  is  outside  the  area.  It  is  possible  that  the 
nearest  point  on  the  boundary  of  the  area  is  not  one 
of  the  corners,  but  is  a  side. 

A  crude  way  of  considering  this  complication  is  to 
examine  each  set  of  sequential  corners,  determine  wheth¬ 
er  the  point  of  closest  approach  of  the  line  to  the 
point  lies  between  the  corners,  and  if  it  does,  and  if 
the  distance  is  less  than  r  .  replace  by  this  dis- 
tance.  mln 


Consider  two  points,  a  and  b.  The  line  between 
these  is  given  by  the  following  parametric  equation: 

Lx  -  ax  +  (bx-ax)  *  Q 
Ly  =  ay  +  (by-ay)  *  Q 


The  distance  between  a  line  point  and  the  grid  point  is 
given  by  the  equation: 

2  2 
dis  (g,L)  =  (gx-Lx)  ,  (gy-Ly)  ) 


=  (gx-ax  +  (bx-ax) Q) 

+  (gy-ay  +  (by-ay) Q)z) 


1/2 
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To  find  the  minimum  point  the  distance  squared  may 
be  differentiated  with  respect  to  Q. 


3dis2 

3Q 


(gx-  (ax+  (bx-ax)  Q) 


+  (gy-ay+ (by-ay) Q) 


(Bx-ax) 

(by-ay) 


=  0  if  and.  only  if 
_  (bx-ax)  .  (qx-ax)  +  (by-ay)  .  (cry-ay) 
(bx-ax) 2  +  (by-ay) 2 


If  Q  >  1  or  Q  £  0  this  line  need  not  be  con¬ 
sidered  since  the  closest  point  is  not  between  a  and 
b.  If  0  <  Q  <  1  then  Q  may  be  substituted  in  equa¬ 
tion  (3)  to  determine  a  minimum  distance  for  that  line. 

By  the  preceding  operations  it  is  possible  to 
determine  both  a  minimum  and  maximum  radius  such  that 
a  ring  centered  at  the  point  (gx,  gy)  would  intersect 
the  vulnerable  area  if  and  only  if  its  radius  r  is 
between  rm^n  and  rmax.  By  examining  each  area  in 
turn,  it  is  possible  to  obtain  the  total  set  of  these 
radii.  These  radii  may  then  be  ordered  and  any  over¬ 
laps  may  be  combined.  Overlaps  occur  when  ftie  maximum 
radius  that  intercepts  one  area  is  larger  than  the 
minimum  radius  that  intercepts  the  next  further  area. 

In  this  case,  a  radii  interval  from  the  minimum  radius 
of  the  closest  area  to  the  maximum  radius  of  the  fur¬ 
ther  area  can  be  obtained. 

Using  these  techniques,  a  set  of  radii  intervals 
for  which  intercept  occurs  may  be  generated  for  each 
grid  point.  These  will  then  be  temporarily  stored 
along  with  the  grid  point  coordinates  on  tape  for  fur¬ 
ther  processing  in  stage  F. 

Stage  E:  Centroid  Distribution  Input.  The  input  in¬ 
formation  about  the  centroid  location  used  up  to  this 
point  in  the  program  was  only  used  to  generate  the 
grid  points.  These  should  have  been  chosen  so  that 
the  data  that  have  been  stored  on  tape  may  be  used  for 
a  large  number  of  specific  centroid  distributions. 

To  further  process  the  data,  a  specific  distribution 
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must  be  input.  A  centroid  distribution  is  determined 
by  its  bias,  its  two  standard  deviations,  and  by  the 
direction  of  the  maximum  standard  deviation  with  re¬ 
spect  to  the  x  axis. 

bcx  =  bias  of  centroid  in  x  direction 
bey  =  bias  of  centroid  in  y  direction 
scu  =  standard  deviation  of  centroid  in 
major  axis  direction 
sev  =  standard  deviation  of  centroid  in 
minor  axis  direction 
ag  =  angle  between  x  axis  and  major 
axis  of  distribution. 

Since  the  test  points  have  X  and  Y  coordinates  it 
is  useful  to  determine  the  variance  along  the  X  axis 
along  the  Y  axis,  and  the  covariance  between  the  two 
axes: 

2  2  2  2 
Vx  =  cos{ag)2  *  sex^  +  sin(ag>2  *  sev^ 

Vy  =  sin(ag;  *  scu  +  cos(ag;  *  sev 

CVyx  =  cos(ag)  .  sin(ag)  *  (scu-scv) 

The  probability  density  at  any  test  grid  point  is 
then  determined  by  the  normal  density  equation: 

den(av,av)=  - ± - *  exp(-(fe*-~bcx)  2  +  (Mr**?)  2 

aeniq..,gyi  2?r  Vx.Vy  expv  V  2*Vx  2*Vy 

+  (qx-bex) (qy-bey) » * 

CVxy  1  '  ’ 

2 

The  area  represented  by  each  grid  square  is  sep  ,  thus 
the  approximate  probability  of  being  within  the  area 
represented  by  a  grid  point  is 

2 

p(gx,gy)  =  sep  *  den  (gx,gy).  (5) 

Stage  F:  Probability  of  Intercept  Given  any  Radius. 

By  using  the  radii  intervals  already  computed,  and  the 
probability  that  a  shot  is  within  any  grid  square,  it 
is  possible  to  generate  the  probability  ‘".hat  an  inter¬ 
cept  occurs  for  any  specified  radius,  Jv-r  further  com 
putation  we  will  produce  a  table  of  probability  of  hit 
for  discrete  radii.  The  table  will  c-.  itain  the  proba¬ 
bility  for  ngd  radii  intervals  with  radii 
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r  =  r  zero  ±{K  .  cut(ngd)  *  sr/ngd) 

ACj  '  =  (cut(ngd)  *  sr/ngd) 

The  information  stored  by  stage  D  can  be  used  to  gen¬ 
erate  the  table  of  Pi(rj),  the  probability  of  inter¬ 
cepting  a  vulnerable  area  with  a  ring  whose  radius  is 
between  r j  and  r j  + Ar  j . 

The  stored  information  contains  gx  and  gy  so 
that  the  probability  of  being  within  the  area  repre¬ 
sented  by  the  grid  point,  P(gx,gy) ,  can  be  calculated 
by  equation  (5).  After  the  Pi(rj)  table  has  been 
initially  zeroed,  the  P(gx,gy)  should  be  added  to  each 
of  the  cells  of  Pi(rj)  where  rj  is  such  that  inter¬ 
cept  occurs.  In  addition,  estimates  will  be  made  of 
the  probability  that  a  grid  point  will  have  interfering 
radii  outside  the  range  of  those  represented  in  the 
matrix,  or  that  the  intercept  point  is  entirely  outside 
the  grid  considered. 

This  stage  produces  the  second  set  of  intermediate 
data,  the  probability  of  intercept  for  each  radius. 

For  design  and  evaluation  of  warheads  these  data  should 
be  output. 

Stage  G:  Intercept  Probability  for  a  Specific  Radii 
Distribution.  The  distribution  of  ring  radii  will  be 
specified  by  an  average  radius,  r  zero,  and  a  standard 
deviation,  sr.  The  probability  of  intercept  may  be  ob¬ 
tained  by  multiplying  the  probability  of  intercept 
given  a  specific  radius  interval  by  the  probability  of 
that  specific  radius  interval  occurring  given  by  the 
normal  function: 


Pr(r.)  = 


•J2tt  sr 

and  summing  over  all  radii.. 


exp  (_(  l£l-r_se,rpj... , , 


2*sr 


As  before,  an  indication  should  be  given  as  to  the 
probability  that  a  radius  is  less  than  the  smallest  or 
greater  than  the  largest  radius  considered. 
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The  answers  generated  from  this  state  are  the  final 
result  of  the  program  and  are  output.  After  this  is 
done,  the  program  attempts  to  read  more  data.  The  data 
may  cause  the  program  to  return  to  Stage  A,  Stage  E,  or 
Stage  G,  depending  on  the  amount  of  new  data  available. 
If  no  new  data  are  read,  the  program  terminates. 

A  feasibility  study  was  made  to  examine  the  utility 
of  this  approach  for  analyzing  complex  warheads  and  com¬ 
plex  vulnerable  regions.  This  study  indicates  that  this 
approach  is  practical,  but  obviously  slower  than  a  tech¬ 
nique  which  leads  to  an  analytic  or  computational  solu¬ 
tion.  Unfortunately,  it  also  indicates  that  the  analy¬ 
sis  of  complex  situations  will  probably  require  the  type 
of  iterative  solution  outlined. 
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